Rotational energy term in the empirical formula for the yrast energies in even-even 

nuclei 
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We show that part of the empirical formula describing the gross features of the measured yrast 
energies of the natural parity even multipole states for even-even nuclei can be related to the 
rotational energy of nuclei. When the first term of the empirical formula, aA -1 , is regarded as the 
rotational energy, we can better understand the results of the previous analyses of the excitation 
energies. We show that the values of the parameters a and 7 newly obtained by considering the 
aA 1 term as the rotational energy of a rigid rotor are remarkably consistent with those values 
extracted from the earlier 'modified' \ 2 analyses, in which we use the logarithms of the excitation 
energies in defining the 'modified' \ 2 values. 
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In a series of works [H, 0, 01, 0, , an empirical formula 
was proposed to represent the gross features of the yrast 
energies E x of the natural parity even multipole states 
including 2 + , 4 + , 6 + , 8 + , and 10 + for even-even nuclei 
throughout the whole periodic table. This formula is 
expressed in terms of the mass number A, the valence 
proton number N p , and the valence neutron number N n 
as 

E x = aA-~< + p exp(-X p N p ) + /3„ cxp(-A„iY„), (1) 

where six free parameters a, 7, j3i, and Aj(i = p, n) are 
fixed so that the experimental excitation energies can be 
fitted by Eq. (fTJ) for each multipole state Q . This empiri- 
cal formula has been quite successful not only in explain- 
ing the main features of the measured excitation energies 
as a function of mass number A but also in reproducing 
the characteristic simple patterns observed in the N p N n - 
plot of the measured excitation energies @, Q • In spite of 
the success of Eq. (Q}, however, it was not clear how this 
simple formula could represent the overall features of the 
yrast energies so well. Very recently, it was suggested in 
Ref . Q that Eq. ([T]) could be approximated as 



E" 



yA 



(2) 



around the doubly mid-shell region, where N p and N n are 
quite large. In this region, the two exponential terms in 
Eq. Q}, P P exp(— XpNp) + /3 n exp(-A„iV„), are very small 
compared with the first term, a A 1 . Ref. Q also sug- 
gested that the two parameters, a and 7, carry the infor- 
mation about the effective moment of inertia. However, 
it was not explained how these parameters could be re- 
lated to the effective moment of inertia. 

In this Brief Report, wc will first show that the first 
term aA~ 7 can be indeed expressed in terms of the ef- 
fective moment of inertia. Then we will redo some of 
the previous analyses in Ref. 3] but with the first term 
aA f fixed as the rotational energy of a rigid rotor, and 
compare the new results with the previous ones. 



It is well known that a nucleus near the doubly mid- 
shell region has a rotational band, which consists of dif- 
ferent total angular momenta J but shares the same in- 
trinsic state. The energy spectrum of the rotational band 
for J 71 —2 + , 4 + , 6 + , • ■ • with the intrinsic angular mo- 
mentum K — can be written as |St| 



E m t{J + ) — 



J(j + i)h 2 
21 



(3) 



where I is the effective moment of inertia of the nucleus. 
Let us assume that a nucleus is a rigid body which has the 
axial symmetry about the intrinsic 3 axis. The moment 
of inertia / r ig of such a rigid body can be expressed as 



/rig = \MRl (1 + i) 



(4) 



where M is the mass of the nucleus given by M = Au, 
u being the atomic mass unit. The distortion parameter 
8 « (i?3 — R±)/Rq is typically 0.2 ~ 0.3 for nuclei with 
150 < A < 188, where R3, R±, and Ro are the radius of 
a nucleus along the intrinsic 3 axis, in the direction per- 
pendicular to it, and the mean radius, Ro = l^^fm, 
respectively 8| . Since the observed moments of inertia are 
smaller than the moment of inertia given by the simple 
form in Eq. ([4]) by roughly a factor of 2 for nuclei around 
the doubly mid-shell region with 150 < A < 1880, we 
may introduce a factor k to take into account the differ- 
ence between the effective moment of inertia / and I Ylg 
of Eq. (@]) so that 

/ = fe/ rig . (5) 
By inserting Eq. |5J into Eq. ([3|) , we immediately get 



E ro t(J 
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(6) 



where 



with 



a' = a J(J+l) and 7' = 5/3 (7) 
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TABLE I: Six parameters a, 7, f3i, and A*(i = p,n) in the empirical formula of Eq. ([T} are listed for three cases. In the upper 
part, four parameters /3j and A; determined by using fixed 7 = 7' = 5/3 and a = a' = aoJ(J + 1) with cto = 65.96 MeV are 
listed: case i). In the middle part, a, f3i, and determined with fixed 7 = 1.40 are listed: case ii). In the lower part, the 
parameters previously extracted are quoted from Table 2 of Ref. [3] : case iii). iVo refers to the number of data points for each 
multipole state. The values of the 'modified' x 2 are a l so listed. 
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a(MeV) 


Qo(MeV) 


/3 P (MeV) 


/?„(MeV) 


Xp 


A n 


x 2 


No 


27 


1.67 


395.76 


65.96 


0.79 


1.09 


0.42 


0.29 


0.157 


557 


47 


1.67 


1319.20 


65.96 


1.12 


1.54 


0.34 


0.24 


0.094 


430 




1.67 


2770.32 


65.96 


1.31 


1.46 


0.32 


0.18 


0.086 


375 


K 


1.67 


4749.12 


65.96 


1.27 


1.34 


0.26 


0.17 


0.060 


309 


10+ 


1.67 


7255.60 


65.96 


1.30 


1.46 


0.23 


0.18 


0.040 


265 


27 


1.40 


89.89 


14.98 


0.82 


1.15 


0.41 


0.28 


0.126 


557 


47 


1.40 


297.87 


14.89 


1.20 


1.67 


0.33 


0.23 


0.071 


430 


67 


1.40 


654.71 


15.59 


1.40 


1.64 


0.31 


0.18 


0.069 


375 


87 


1.40 


1155.90 


16.05 


1.34 


1.50 


0.26 


0.15 


0.053 


309 


10+ 


1.40 


1702.79 


15.48 


1.34 


1.64 


0.22 


0.15 


0.034 


265 


27 


1.34 


68.37 


11.40 


0.83 


1.17 


0.42 


0.28 


0.126 


557 


47 


1.38 


268.04 


13.40 


1.21 


1.68 


0.33 


0.23 


0.071 


430 


67 


1.38 


598.17 


14.24 


1.40 


1.64 


0.31 


0.18 


0.069 


375 


87 


1.45 


1438.59 


19.98 


1.34 


1.50 


0.26 


0.15 


0.053 


309 


10+ 


1.47 


2316.85 


21.06 


1.36 


1.65 


0.21 


0.14 


0.034 


265 



Eq. ([6]) shows that the first term aA^ 1 of the empir- 
ical formula in Eq. ([1]) can be derived and identified as 
the rotational energy of a deformed nucleus. However, 
this interpretation of the first term aA~ J is valid only 
for nuclei near the doubly mid-shell region, but not for 
closed shell nuclei or for light nuclei with A < 30, most 
of which do not have rotational bands. For closed shell 
nuclei, on the other hand, the first term aA^ 1 is negligi- 
ble compared to the other two exponential terms, which 
will be shown later in Fig. 3. Thus if we accept this inter- 
pretation to this extent, Eq. ([6]) shows that 7 in Eq. ([1]) 
can be identified as a constant 7' close to 1.67 for all of 
the even multipole states. Note that the values of 7 in 
Eq. {1} extracted earlier and shown in Table 2 of Ref. 3] 
range from 1.34 to 1.47. These previous values of 7 are 
listed again in the lower part of Table |U The average 
of these previous 7 values is 1.40 which is remarkably 
close to 7' = 5/3, in spite of the fact that the values of 
7 were extracted earlier without any constraints on 7 or 
consideration presented in Eqs. (O^®. 

Let us repeat the x 2 analyses as was done in Ref. [3| but 
with a constant average value of 7 = 1.40. In the fitting 
procedure we use as in Ref. 0] the logarithm of the ex- 
citation energies E%. (i) and E^ xp (i) since the excitation 
energies are spread over a wide range. This 'modified' % 2 
analysis is different from the conventional x 2 analysis in 
that the small values of the excitation energies are em- 
phasized. By defining R E (i) = log [E!f '(i)] -log [E° x P(i)], 
the five parameters a, Pi and Ai are fixed to minimize 
X 2 = jf- J2^=i I^e(*)] 2 , where N is the total number of 



the data points considered for the corresponding mul- 
tipole state, a, Pi, and Ai extracted in this way are 
given in the middle part of Table Q] The values of 
a and ao = a/J(J + 1) determined with 7 = 1.40 
are somewhat different from those of a and ao in the 
lower part of Table HI On the other hand, the values 
of Pi and Ai in the middle part of Table Q] are still very 
close to those values in the lower part for all the multi- 
pole states. This shows that the two exponential terms 
P p exp(— XpN p ) + p n exp(—X n N n ) remain rather robust 
even if a and 7 change and thus that the two exponential 
terms can be well separated from the rotational energy 
term. 

Let us now fix a and 7 as given by Eq. ([7]) for all J and 
redo the 'modified' x 2 analyses for 2+ to 10+ states. In 
the analyses we use k = 1/2 and the distortion parameter 
S = 0.3 determined for 174 Yb (N p = 12, N n = 22) which 
is a nucleus in the doubly mid-shell region Q . Then the 
value of ao becomes 65.96 MeV. In TableUwe show that 
the 'modified' x 2 values obtained by using these fixed val- 
ues of a' and 7' of Eqj[7]) turn out to be larger than those 
previously obtained H by 25% (for 2+), 32% (for 4+), 
25% (for 6+), 13% (for 8+), and 18% (for 10+). Never- 
theless, the excitation energies calculated by using fixed 
ao = 65.96 MeV and 7' = 5/3 appear almost identical to 
those obtained earlier in 0]. 

To compare the different values of a listed in Table QJ 
we plot in Fig. Q]the values of a and ao against J for three 
cases: i) 7 = 7' = 1.67 (circles), ii) 7 = 1.40 (squares), 
and iii) 7 from Ref. 3] (triangles), a's for all three cases 
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FIG. 2: (Color online) The first term of Eq. |TJ for the 2+ 
states in even-even nuclei are plotted against the mass number 
A. The solid curve represents the values of qA~ 7 calculated 
with a and 7 from Ref. 0]. The dashed curve denotes the 
values of a.A~"' with 7 = 1.40 and a determined by the \ 2 
fitting. The dotted curve denotes aA ' when ao and 7 are 
fixed as 65.96 MeV and 1.67, respectively. 



FIG. 1: (Color online) The values of a and ao in Table U are 
plotted for the following three case: i) 7 = 1.67 (circles), ii) 
7 = 1.40 (squares), and iii) 7 from Ref.[3| (triangles). 



in Fig. QJa) change as a quadratic function of J. Thus 
we plot in Fig. [TJb) a = a/J{J + 1). Figure [TJb) 
shows that ao's are indeed constant for all three cases. 
It is remarkable to see that the values of ao (triangles) 
obtained from those of a extracted previously in Ref. 3] 
turn out to be also almost constant. Figure Q] shows that 
the values of a for the cases ii) and iii) are very close 
to each other but are smaller than those for the case i) 
roughly by a factor 4. This factor is to compensate for 
the changes in aA~ n due to the changes in 7 from about 
1.40 to 1.67. Note that the ratio A~ 1A0 /A~ 1S7 = A - 27 
is about 4 for 17i Yb. 

In Fig.[2]the value of aA~~* term with a and 7 as given 
in Ref. [3j is plotted by the solid curve, and that with 
7 = 1.40 is plotted by the dashed curve. They agree very 
well. The values of E mt (J + ) of Eq. (O with a' and 7' 
given in Eq. are plotted by the dotted curve, which 
also agree with the other curves quite well for larger 
A. This shows that the values of a compensate for the 
changes in the values of 7 so that aA~ 7 remains more 
or less the same. Note, however, that the dotted curve 
deviates from the other two curves for small values of A. 
It is because we use the logarithm of energies rather than 
the energies themselves for the 'modified' \ 2 calculations. 
From the x 2 analyses with ao = 65.96 MeV and 7 = 1.67 
we obtain new parameters f3i and Xi listed in the upper 



part of Table HI The new parameters of /3j and Xi are 
still very close to those values in the middle and lower 
parts. As mentioned earlier, the parameters of (3i and Xi 
are almost the same for three cases of i) ~ iii), implying 
that the two exponential terms are rather independent of 
the rotational energy term a A 1 . 

In Fig.[3] we show the yrast energies for each state of 
2+ (a), 4+(b), 6+(c), 8+(d), and 10+ (e) against the mass 
number A. For each J v state the experimental yrast en- 
ergies are plotted in the upper part and are compared 
with the calculated values in the lower part obtained by 
using the parameters given in the upper part of Table 
HI The measured excitation energies of the first 2+ are 
quoted from the compilation of Raman et al. 0] and those 
of the first 4+ ~ 10+ are extracted from Ref. 0. The 
data points are connected by the solid lines along the 
isotopic chains. Figure[3] shows that Eq. ([TJ) with fixed a' 
and 7' in Eq. ([7]) can reproduce the overall trends of the 
measured yrast energies for all of the states. To show the 
discrepancies between the experimental energies and the 
calculated energies, we plot the ratios £"^ al /i?° xp against 
the mass number A in Fig.UJ (We plot the ratios only 
for the 2+ and 8+ states, because the ratios for the other 
states are quite similar.) The ratios are close to unity for 
nuclei near the doubly mid-shell region. For nuclei in the 
closed shell region where our interpretation of aA 1 is 
not valid, the ratios can be as big as 2. For light nuclei 
to which our interpretation is not applicable E^ 1 is often 
much larger than E° xp . This is because the 'modified' 
X 2 analyses put more emphasis on the smaller values of 
the energies and thus a A" 1 with 7 = 1.67 (dotted curve 
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FIG. 3: (Color online) The measured excitation energies of the 2j*", 4J~ , 6^, 8^, and 10^ states in even-even nuclei are compared 
with our results. The measured data are connected by the solid lines along the isotopic chains. The upper parts of (a)-(e) show 
the measured energies 0, [nj, while the lower parts show the calculated energies by using the parameter sets in the upper part 
of Table U The solid curve and the circles in the bottom part of (a) show the contribution to the calculated energies of the 2^~ 
states from the first term and two exponential terms. 



in Fig. [5]) overshoots aA~^ with 7 = 1.40 (dashed curve 
in Fig. [2]) for small values of A. Finally we show in the 
bottom part of Fig. 3 (a) the separation of the calculated 
energies into the first term aA~ 7 and the two exponen- 
tial terms. The first term aA~^ (the solid curve) is the 
major contribution in the yrast energies in the doubly 
mid-shell region while it is negligible compared to the 
two exponential terms (circles) in the closed shell region. 

In summary, we have shown that the term aA~ 7 can 
be obtained by considering the moment of inertia of a 
deformed nucleus. The yrast energies calculated with 
constant a' and 7' can describe the main features of the 
data. It is remarkable that the values of aA ' extracted 
earlier in Ref. Q agree well with the values of a'A~ 7 ob- 
tained from the rotor model and that the previous values 
of a divided by J(J+ 1) are almost constant as expected 



from the rotor model. It is also seen that the parame- 
ters Pi and Xi newly extracted with constant a' and 7' 
in this work are very consistent with those parameters 
previously obtained in Ref. 3]. It shows that the em- 
pirical formula can be well separated into the rotational 
energy term aA 1 and the two exponential terms which 
are thought to be related to the shell effect. 
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FIG. 4: The ratios of the calculated excitation energies to the 
measured energies for the 2± and 8f , E^ 1 / ' E^ v , are plotted 
against the mass number A. 
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